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ABSTRACT 

Inspired by increase of population of 7-ray emitting pulsars by the Fermi 
telescope, we perform a population study for 7-ray emitting canonical pulsars. 
We use a Monte-Carlo technique to simulate the Galactic population of neutron 
stars and the radio pulsars. For each simulated neutron star, we consider the 7- 
ray emission from the outer gap accelerator in the magnetosphere. In our outer 
gap model, we apply the gap closure mechanism proposed by Takata et al., in 
which both photon-photon pair-creation and magnetic pair-creation processes are 
considered. Simulating the sensitivities of previous major radio surveys, our sim- 
ulation predicts that there are ~ 18 — 23 radio loud and ~ 26 — 34 7-ray-selected 
7-ray pulsars, which can be detected with a 7-ray flux F 7 > 10~ 10 erg/cm 2 s. Ap- 
plying the sensitivity of the six-month observation of the Fermi telescope, 40-61 
radio-selected and 36-75 7-ray selected pulsars are detected within our simulation. 
We show that the distributions of various pulsar parameters for the simulated 7- 
ray pulsars can be consistent with the observed distribution of the 7-ray pulsars 
detected by the Fermi telescope. We also predict that ~ 64 radio-loud and ~ 340 
7-ray-selected pulsars irradiate the Earth with a flux F 7 > 10 -11 erg/cm 2 s, and 
most of those 7-ray pulsars are distributing with a distance more than 1 kpc and 
a flux F 7 ~ 10" 11 erg/cm 2 s. The ration between the radio-selected and 7-ray- 
selected pulsars depend on the sensitivity of the radio surveys. We also discuss 
the Galactic distribution of the unidentified Fermi sources and the canonical 
7-ray pulsars. 

Subject headings: pulsars:general-gamma rays:theory-radiation mechanisms :non- 
thermal-stars:neutron 
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1. Introduction 

The Fermi 7-ray telescope has drastically increased population of 7-ray emitting 
pulsars after its launch in June, 2008. Abdo et al. (2010a) reported the first Fermi pulsar 
catalog, including 21 radio selected pulsars, 16 7-ray-selected pulsars, the Geminga and 
8 millisecond pulsars (see also Abdo et al. 2009a,b), while the pulsed radio emissions are 
reported from several 7-ray-selected pulsars (Camilo et. al. 2009). Saz Parkinson et al. 
(2010) have also reported new eight gamma-ray pulsars discovered in blind searches of 
Fermi- LAT data. These Fermi observations for the 7-ray emitting pulsars have started to 
break the bottleneck of study of the high-energy radiation from the pulsar magnetosphere. 

The particle acceleration and resultant high-energy emissions from the pulsar 
magnetosphere have been studied with polar cap model (Ruderman & Sutherland 1975; 
Daugherty & Harding 1982, 1996), slot gap model (Arons 1983; Muslimov & Harding 2003; 
Harding et al. 2008) and outer gap model (Cheng, Ho & Ruderman 1986a,b; Hirotani 2008; 
Takata, Wang & Cheng 2010; Wang, Takata & Cheng 2010). All emission models assume 
that the electrons and/or positrons are accelerated by the electric field along the magnetic 
field lines and that the accelerated particles emit the 7-rays either the curvature radiation 
or the inverse- Compton processes. The polar cap model assumes the emission region close 
to the stellar surface and above the polar cap region, while the outer gap and slot gap 
models assume an acceleration region extending to the outer magnetosphere. 

Properties of the pulse profiles and spectra of the 7-ray emission measured by the 
Fermi can be used to discriminate between the different models. For example, most of 
the observed pulse profiles of the 7-ray pulsars detected by Fermi show the double peak 
structure with a phase-separation of ~ 0.4 — 0.5, indicating a high-altitude emission 
from the outer gap or slot gap region is favored more than the low altitude emission 
from the polar cap accelerator (Romani & Watter 2010; Venter, Harding & Guillemot 
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2009). Furthermore, in the first Fermi pulsar catalog (Abdo et al. 2010a), the spectral 
fits have been done assuming exponential cut-off function. It is found that the observed 
spectrum for the Vela pulsar (Abdo et al. 2009c, 2010b) is well fitted by a power law plus 
hyper-exponential cutoff spectral model of the form dN/dE oc E~ a exp[— (E/E c ) b ] with 
b < 1 and E c ~ 1.5 GeV. This implies that an emission process in the outer magnetosphere 
is favored more than that near the stellar surface, which predicts a super exponential cut-off 
(6 > 1) due to the magnetic pair-creation prosses of the 7-rays. In addition, the detection of 
the radiation above 25 GeV bands associated with the Crab pulsar by the MAGIC telescope 
has also implied the emission process in the outer magnetosphere (Aliu et al. 2008). 

The increase of population of the 7-ray emitting pulsars also allow us to perform a 
more detail study for the high-energy emissions from the pulsars. Wang et al. (2010) fit the 
observed phase-averaged spectra using a two-layer outer gap model, and then investigated 
the relation between the 7-ray luminosity (L 7 ) and the spin down power (L S( i) for 39 7-ray 
emitting pulsars including 8 millisecond pulsars. They suggested that the relation can be 
expressed as L 7 oc L^ d with (3 ~ for L s d > 10 36 erg/s, while (3 ~ 0.5 for L s d < 10 36 erg/s. 
They further argued that the relation between the 7-ray luminosity and the spin-down 
power implies that the outer gap is closed by the photon-photon pair-creation process 
for L s d > 10 36 erg/s (Zhang & Cheng 1997), while by magnetic pair-creation process for 
L s d < 10 36 erg/s (Takata et al. 2010). 

For the outer gap model, the population studies were done to compare the model 
prediction with the results of EGRET instrument (Cheng & Zhang 1998; Zhang et al. 
2004; Cheng et al. 2004). Applying the photon-photon pair-creation model for gap closure 
process, they found that the distributions of the various pulsar parameters (e.g. rotation 
period, spin down age) were consistent with the observations. Gonthier et al. (2002) and 
Story, Gonthier and Harding (2007) studied the population of 7-ray pulsars with the polar 
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cap and slog gap accelerator models to predict the results using the Fermi observations. 
Furthermore, Grenier & Harding (2010) has studied the distribution of the 7-ray pulsars 
predicted by the outer gap model of Zhang et al. (2004) and for the slot gap model of 
Muslimov & Harding (2003; 2004) to compare with the Fermi results. They have argued 
that both outer gap and slot gap model may not explain all features of the Fermi pulsars; 
for example, both models predict too few young pulsars compared with the Fermi pulsars. 

In this paper, we will further develop the population study for the 7-ray emitting 
pulsars to compare with the results of six-month Fermi observations. In particular, we 
will consider the 7-ray emissions of the canonical 7-ray pulsars with the outer gap model. 
We will apply the gap closure model proposed by Takata et al. (2010), in which the 
photon-photon pair-creation process and the magnetic pair-creation process are taken into 
account. As usual, we will perform a Monte-Carlo simulation for the Galactic population 
of the neutron star (section [2]). In section [2~5l we will discuss our 7-ray emission model 
and the thickness of the outer gap with the pair-creation processes. In section [3], we 
will compare the simulated number of the 7-ray pulsars detected on the Earth and its 
distributions with the Fermi observations (section 3.1 and 3.2), and we will predict the 
future observations (section 3.3). In section HI we will also discuss ratio of the 7-ray-selected 
and radio-selected 7-ray pulsars with different radio surveys (section 4.1), and will compare 
Galactic distributions of the simulated 7-ray pulsars and of the Fermi unidentified sources 
(section 4.2). In section [51 we will summarize the results of our population study. 

2. Theoretical model and Monte Carlo Simulation 

In this section, we will describe our Monte Carlo simulation and 7-ray emission model 
for the canonical pulsars. First we perform the Monte Cairo simulation for the Galactic 
neutron stars and for the radio emissions from the pulsars, which have been developed by 
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several authors (e.g. Sturner & Dermer 1996; Cheng & Zhang 1998, Gonthier et al. 2002). 
With this simulation, we will be able to obtain a simulated population of the radio pulsars 
consistent with the observations. Then we will consider the 7-ray emissions from the outer 
gap accelerator of each simulated neutron star. 

We note that we will sample and discuss the 7-ray emissions from the simulated 
neutron stars with a characteristic age younger (r = P/2P) than 10 Myr. This is because 
the present known 7-ray emitting canonical pulsars have their characteristic spin down age 
less than 5 Myr. Also, we will find that the present 7-ray emission model and the simulation 
predict that the canonical pulsars with age older than 10 Myr are hard to produce the 
observable 7-ray emissions in their magnetosphere. 

2.1. Period and magnetic field 

We assume that the neutron stars are born at a rate of 1~2 per century. The initial 
period of the new born neutron star is not well predicted by theory and observation, except 
for several cases. For the Crab pulsar, the initial period is estimated at P$ ~ 19 ms. Such 
short birth rotation period is also expected by the existence of young 16 ms pulsar PSR 
J0537-06910 (Marshall et al. 1998). Although most of the new pulsars may have initial 
period Pq ~ 20 ms, there is good evidence that some pulsars were born with a longer initial 
period (e.g. P ~ 62 ms of PSR J1811-1925, Kaspi et. al., 2001). However, we found that 
the distribution of initial rotation period does not affect much to the following results of the 
simulation (e.g. the distribution of the radio pulsars and of 7-ray pulsars). For example, 
we compared the results of two simulations with different distribution of the initial period; 
one is that all simulated pulsars are randomly distributed in the range Pq = 20 — 30 ms, 
and other is that 70 % of pulsars are distributed at P = 20 — 30 ms and 30 % are 
Po = 20 — 100 ms. As we will describe in later, because the pulsars younger than 5 My can 
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become the 7-ray emitting pulsars, we will sample the pulsars younger than 10 My. For 
the population of the relatively younger pulsars, we could see that the populations of two 
simulations are very similar to each other. The difference becomes significant if we assume 
> 50 % are randomly distributed in P = 20 — 100 ms. In this paper, therefore, because 
we do not know well exact distribution of initial rotation period, we randomly choose the 
initial period in the narrow range Pq = 20 — 30 ms, such like the assumption in the previous 
Monte-Carlo studies. 

We assume a Gaussian distribution in log 10 B s for the initial strength of the magnetic 
field on the stellar surface, 



In this study we apply log 10 Bo = 12.6 and ob = 0.1, which produce a consistent distribution 
of the magnetic field of the known radio pulsars derived from the dipole radiation model 
(Manchester et al. 2005). Because we simulate the neutron stars with a spin down age 
younger than r < 10 Myrs, we do not consider the evolution of the magnetic field of each 
simulated neutron star. The evolution of the stellar magnetic field may be important for 
the neutron star with an age larger than r > 10 Myrs (Goldreich & Reisenegger 1992; 
Hoyos, Reisenegger & Valdivia 2008). 

With a constant stellar magnetic field with time, the period evolves as 



where we assume the pure dipole radiation, R s is the stellar radius and / is the neutron 
star momentum of inertia. We apply R s = 10 6 cm and / = 10 45 g • cm 2 in this study. The 
time derivative of the rotation period is calculated from 




(1) 




(2) 
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2.2. Initial spatial and velocity distributions 

Following the studies done by Paczynski (1990) and Sturner & Dermer (1996), we 
assume the birth location of the neutron star is described by the following distributions, 



Pr{R) 



R 2 



p z (z) = —e^'^, (4) 

where R is the axial distance from the axis through the Galactic center perpendicular to 
the Galactic disk and z is the distance from the Galactic disk. In addition, a,R = 20 kpc, 
i? cxp = 4.5 kpc, or = [1 - e -- Bm «*/ fl «*(l + Rmax/Rexp)}" 1 with R ma x = 20 kpc and 

z exp = 75 pC. 

For the distribution of the initial velocity of simulated neutron star, we refer the 
result of the study for the radio pulsars done by Hobbs et al. (2005). They found that 
the distribution of the three-dimensional velocity of their sample is well described by a 
Maxwellian distribution with a characteristic width of a v = 265 km/s, namely, 

PM = yj\ (5) 

We apply this form for the distribution of the initial kick velocity due to the super nova 
explosion. In this study, we randomly choose the direction of this kick velocity of the 
neutron star in the three-dimensional space. For the azimuthal component of the velocity, 
we add the circular velocity due to the Galactic gravitational potential field at the birth 
position of the neutron star. The circular velocity due to the Galactic potential is calculated 
from 



■ d<$> sph d<$> dis d$ h 



- 1/2 

OR ' OR ' OR 1 
where § sp h, &dis and $^ are spheroidal, disk and halo components of the Galactic 
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gravitational potential, which are given by Pacynski (1990) as 

i / \ GMi 



(7) 



where a sp h = 0, b sp h = 0.277 kpc and M sp h = 1.12 x 1O 1O M for the spheroidal component 
a>dis = 3.7 kpc, bdis = 0.20 kpc, and = 8.07 x 1O 1O M for the disk component, and 



GM C 



In 1 + 



' r 1 

+ — tan" 1 



where r c = 6.0 kpc M c = 5.0 x 1O 1O M for the halo component 



2.3. Equation of motion 

To obtain current position of each simulated neutron star, we solve the equation of 
motion from its birth to the current time. The equation of motion is given by 



dR 2 vl 9$ 



dt 2 R dR 



(9) 



dz 2 <9$ 



(10) 



dt 2 dz ' 
and 

RV0 = constant, (11) 

where $ tot = § sp h + &dis + and is the azimuthal component of the velocity. The 
Lagrangian in units of energy per unit mass is described by 

L= -W (12) 

which is used for maintaining an accuracy of the integration of the trajectory described by 
equations (Q and fTTOj) . We maintain the accuracy of one part per 10 5 during the trajectory 
from its birth to the current time. 
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2.4. Radio emission and detection 

Because detail process of the radio emission in the pulsar magnetosphere has not been 
understood well, we apply an empirical relation among the radio luminosity, the rotation 
period and the period time derivative (Cheng & Zhang 1998; Gonthier et al. 2002). We 
assume that the radio luminosity at 400 MHz, which is defined by L 40 o = d 2 S/too with 
d being the distance and S400 observed flux, for each simulated neutron star follows the 
distribution (Narayan & Ostriker 1990) 

p Li00 = 0.5AV, (13) 

where A = 3.6[log 10 (i:4oo/ < L 400 >) + 1.8] with log < L 400 >= 6.64 + |logi (P/P 3 ), 
and L400 is the luminosity in units of mJy kpc 2 . The radio flux on the Earth is given by 
S i00 = L4 00 /cP. We scale the simulated 400 MHz luminosity to the observation frequency of 
each survey using a typical photon index -2. 

The radio emissions from the pulsars may not point toward the Earth. We apply the 
radio beaming fraction described as (Emmering & Chevalier 1989) 

f r {uj) = (1 — cosw) + (7r/2 — oj) + sinu;, (14) 

where w is the half-angle (in radian) of the radio emission cone and < f r < 1. We apply 
the half-angle of the radio emission cone studied by Kijak & Gil (1998, 2003), 

u KG ~ 0.02r^P~ 1 / 2 , (15) 

where 

r _ An, .-0.26 p0.07 p0.3 

where -P-15 is the period time derivative in units of 10 -15 and vghz is the radio frequency 
in units of GHz. In the Monte-Carlo simulation, the fraction f r of the pulsars with same 
width uj can emit the radio beam toward the Earth. 
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The minimum detectable radio flux (that is, sensitivity) of a particular radio survey is 
usually express as (Edwards et al. 2001), 



r, _ C t hres[T r ec + T s k y (l,b)] W . . 

where C t hres is the receiver detection threshold S/N, T rec is the receiver temperature, 
T s k y (l,b) is the sky temperature with / being Galactic longitude and b Galactic latitude, 
G is the telescope gain, Bbd is the band width, is the integration time and W is the 
observed pulse width. The sky temperature at the frequency v is given by (Johnston et al. 
1992) 

nn , S or f 2 ^ ] /408 MHz\ 2/6 „ 

The observed pulse width is given by 



W 2 = Wl + T 2 samp + T 2 DM + T 2 cat , (18) 

where Wo ~ 0.05P is the intrinsic pulse width, r samp is the receiver sampling time, tdm is 
the frequency dispersion of the radio pulse as it passes the interstellar medium, and r scat is 
the frequency dispersion due to the interstellar scattering. The dispersion t^m has the form 

r M/ = 8.3xlO«DM(^)(^-)- 3 ms (19) 

where 5v is the frequency channel bandwidth of the receiver and DM (in units of pc cm -3 ) 
is the dispersion measure, which is calculated from DM = J n e (R, Z)ds, where s is the 
distance to the source and n e (R, Z) is electron number density. The Galactic distribution 
of the electron number density is calculated from the model investigated by Cordes & Lazio 
(2002). The dispersion due to the interstellar scattering is obtained by 

/400 MHz\ 4 ' 4 

Tscat{v) = r scatA 00 I I HIS (20) 



where r scatim = l0- 4 - 61 + L14I °s" DM + iQ-^+^eiog^ dm^ 
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In Table 1, we list the radio surveys that we applied in the present simulation and its 
receiver parameters in equations ffTE^ -f liZUj) . Also, Figure [T] plots the sensitivity described by 
equation (fTT5j) for each pulsar survey as a function of the period. In Figure [TJ we used the sky 
temperature of T s f. y = 150 K at 408 MHz and the dispersion measure of DM=200 cm _3 pc. 
We note that although new radio pulsars have been reported with the search of the pulsed 
emissions for the specific sources (e.g. Camilo et al 2009 for Fermi sources), the majority 
of the present radio pulsars have been discovered by the surveys listed in Table 1. Therefore 
we do not take into account such pulse search for the individual sources. 



In this section, we describe our 7-ray emission model applied for each simulated pulsar. 
We consider the 7-ray emission from the outer gap accelerator (Cheng, Ho & Ruderman 
1986a,b; Zhang & Cheng 1997; Takata et al. 2010). In the outer gap, the charged particles 
(electrons and positrons) are accelerated up to a Lorentz factor of T > 10 7 by the electric 
field along the magnetic field line. The accelerated particles can emit several GeV 7-ray 
photons through the curvature radiation process. Assuming the force balance between the 
electric force and the radiation drag force, the Lorentz factor is described by 



where R c is the curvature radius of the magnetic field lines and En is the accelerating 
electric field. The accelerating electric field in the outer gap is approximately described by 
(Cheng et al. 1986a,b; Cheng, Ruderman & Zhang 2000) 



where Ri c = cP/2ir is the light radius and / gQ p( r ) is so called fractional gap thickness in the 
poloidal plane. The power of the curvature radiation emitted by the individual particle is 



2.5. 7-ray emission model 




(21) 



E \\ ( r ) 



B(r)f gap (r) 2 R lc 

R c 



(22) 
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written as 

P c (E 1 ,r) = ^-F(x), (23) 
where x = E 1 /E c with E c = 3hcT 3 /4:7rR c and 

POO 

F{x) =x K s/3 (t)dt, 

J X 

where -K5/3 is the modified Bessel function of the order 5/3. If the 7-ray beam points toward 
an observer, the observer will measure the phase-averaged spectrum of (Hirotani 2008), 

^~N(r,)R c (r,)P c (E,,r^, (24) 

where N is the typical particle number density, represents the typical radius at emission 
positions from which the emissions can be measured by the observer and d is the distance 
to the source. In addition, 8 A is the cross section of the gap perpendicular to the magnetic 
field lines and is estimated from 8 A ~ [B s /B(r^)]Rp898(f), where R p is the polar cap radius, 
and 89 and 8(f> is gap width in the co latitude and azimuthal directions measured on the 
stellar surface, respectively. For the typical case of the outer gap model, we apply the 
Goldreich- Julian value, N = B/Pce, for the number density, and the curvature radius of 
R c = Ri c /2. We assume that the outer gap is typically extending to the azimuthal direction 
with 8(f) ~ 71. The integrated energy flux between 100 MeV and 300GeV can be calculated 
from 

/•300GeV 7p 

-^7,>100 MeV — / 

JlOOMeV "-^7 

In the paper, we assume that the gap current is order of Goldreich- Julian value over the 
full width as zeroth order approximation. The detailed calculation of the outer gap model 
(e.g. Takata & Chang 2007) shows that the pair-creation process produces the distribution 
of the current density in the direction perpendicular to the magnetic field lines. On the 
other hand, recent two-layer outer gap model (Wang et al. 2010) predicts that total current 
in the gap is order of the Goldreich- Julian value. Therefore we expect that as long as we 
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discuss the total power of the 7-ray radiation from the outer gap accelerator, the uniform 
current distribution with the Goldreich- Julian value would not be bad approximation. The 
distribution of the current in the gap will be more important to discuss the spectral shape 
above 100 MeV bands. 

The beaming fraction / 7 , which is defined by ratio of the solid angle of the 7-ray beam 
to 47r, affects number of the detected 7-ray pulsars. However, it is complicated to describe 
the beaming fraction with a simple function. This is because the beaming fraction should 
be different for different inclination angle and for different pulsars. For the outer gap model, 
however, the most of emissions concentrate within the direction 90° ± 35° measured from the 
rotation axis, as we can see in the photon mappings (e.g. Cheng, Ruderman & Zhang 2000; 
Takata & Chang 2007), indicating the beaming fraction is roughly / 7 ~ 35°/90° ~ 0.4. 

Let's discuss the fractional gap thickness f gap , which determines the 7-ray flux of 
equation (j24"|) . Recently, Wang et al. (2010) studied the relation between the 7-ray 
luminosity and the spin down luminosity for the 7-ray pulsars detected by the Fermi 
telescope. They found that the 7-ray luminosity can be described by 



They further argued that because the typical 7-ray luminosity is described by L 7 ~ fg ap L sd , 
the change of the relation between L 1 and L sd indicates the gap closure mechanism in the 
poloidal plane switches at L sd ~ 10 36 erg/s. 

The relation that L 7 oc L sd with (3 ~ in equation (1261) may be understood with 
the gap closure mechanism proposed by Zhang & Cheng (1997). They argued that the 
photon-photon pair-creation process between the 7-rays emitted in the outer gap and the 
X-rays from the stellar surface controls the gap activities. They discussed that the X-rays 
from the heated polar cap by the return particles will collide with 7-rays in the outer gap 




/3 ~ 0, for L sd > 10 36 erg/s 
P ~ 0.5, for L sd < 10 36 erg/s. 



(26) 
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and provide a self-consistent mechanism to restrict the fractional size of the gap. They 
estimated the typical gap size near the light cylinder as 

f zc = 5.5P 26 / 21 i?r 2 4/7 , (27) 

where B u is the polar magnetic field strength in units of 10 12 G. With this model, the 7-ray 
luminosity and the spin down power is related as 

L 7 ~ 1.3 x \^B^L% erg/s, (28) 

where we used L s d,m is the spin down power in units of 10 36 erg/s. 

On the other hand, the relation that L 7 oc L^ d with (3 ~ 0.5 in equation (1261) may be 
explained by the gap closure model proposed by Takata et al. (2010), in which the outer gap 
is closed by the magnetic pair-creation process near the stellar surface. They argued that 
the returning particles, which were produced by the photon-photon pair-creation process 
in the gap, will emit the photons with an energy m e c 2 /otf ~ lOOMeV by the curvature 
radiation near the stellar surface. They discussed that those lOOMeV photons can become 
pairs by magnetic pair creation process and that these secondary pairs can continue to 
radiate several MeV photons via the synchrotron radiation. The photon multiplicity is 
easily over 10 5 per each incoming particle. For a simple dipole field structure, all pairs 
should move inward and cannot affect to the outer gap accelerator. However they argued 
that the existence of strong surface local field (e.g. Ruderman 1991, Arons 1993) has 
been widely suggested. In particular if the field lines near the surface, instead of nearly 
perpendicular to the surface, are bending sideward due to the strong local field, the pairs 
created in these local magnetic field lines can have an angle bigger than 90°, which results 
in an outgoing flow of pairs. In fact it only needs a very tiny fraction (1-10) out of 10 5 
photons creating pairs in these field lines, which are sufficient to provide screening in the 
outer gap when they migrate to the outer magnetosphere. With this model, they estimated 
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the fractional gap thickness when this situation occurs as 

f m = 0.8KP 1 / 2 , (29) 

where K ~ B^ n 2 12 s 7 is the parameter characterizing the local parameters, B m V2 . and Sj are 
the local magnetic field in units of 10 12 G and the local curvature radius in units of 10 7 cm, 
respectively. They argued that K ~ 2 for the canonical pulsars, while K ~ 15 for the 
millisecond pulsars. This gap closure model predicts that the 7-ray luminosity is related 
with the spin down power as 

L 7 ~ 1.1 x 10™K*B 3 £l% erg/s. (30) 

We can see that the fractional gap thickness fzc is smaller than f m for the younger 
pulsar, implying the gap thickness is determined by the photon-photon pair-creation 
process, while f m < f zc for the older pulsars, implying the magnetic pair-creation process 
may control the gap thickness. Equating fzc and f m , we obtain the critical spin down 
power that 

L sd , c ~ 10 36 K- 168 / 31 i?r 2 34/31 erg/s, (31) 

which may explain the switching positions in equation ( 126]) argued by Wang et al. (2010). 
In this paper, we will investigate the population predicted by the 7-ray emission model 
with the switching of the gap closure mechanism at L s d )C , because the population with the 
emission model without the switching (that is f gap = f zc ) have been investigated in the 
previous studies (e.g. Chang and Zhang 1998; Zhang, Zhang and Cheng 2000; Zhang et al. 
2004). 

Finally we would like to discuss the maximum fractional thickness, f ma x, f° r the active 
outer gap accelerator. Zhang & Chang (1997) argued that the pulsar with the fractional 
gap thickness larger than unity, fzc > 1, is not active, because the pairs are not created in 
the gap by the photon-photon pair-creation process. However, it may be possible that the 
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maximum gap thickness is smaller than unity. For example, outer gap accelerator can exist 
between the last-open field lines and the critical field lines, on which the Goldreich- Julian 
charge density is equal to zero at the light cylinder. Therefore, we may define the maximum 
gap thickness as 

fcrit = ^ j (32) 



9 p 

where 9 P and 9 C are polar angles of the last-closed filed line and of the critical field line on 
the stellar surface, respectively. For the pure dipole field, we obtain 



9v 



a + sin" 



sm{U k 



a 



— sin di c 
J~£ic 



1/2 



and 



9 C = a + sin 1 



sm(U„ — a) 



D \ 1/2 

— sin 9, 

tile 



(33) 



(34) 



respectively, where a is the inclination angle, 



tan 



-3 - V9 + 8 tan 2 



a 



4 tan a 



and 



9, 



tan 



_! / 3 tan a + V8 + 9 tan 2 



a 



In the realistic situation, we expect that the polar cap shape and the magnetic field 
configuration in the pulsar magnetosphere are more complicated and the simple dipole 
formulae for 9 p and 9 C may not be applied. In this paper, therefore, we examine two extreme 
cases that f max =l and f max = (9 P — 9 C )/9 P with the dipole field, where we randomly choose 
the inclination angle a. These results of the two case may give a range of uncertainty of the 
present theoretical predictions. 



3. Results 



In this paper, we sample the pulsars with a characteristic age smaller than r < 10 Myr, 
because (1) we are interesting in the population study for the 7-ray emitting canonical 
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pulsars and (2) the present canonical 7-ray pulsars have been detected with a characteristic 
age younger than 5 Myr. In fact, although we run the simulation to obtain more than 200 
times of the present radio pulsars, we could not find any detectable canonical 7-ray pulsars 
with a characteristic age older than 10 Myr. 

We simulate 3 x 10 6 of neutron stars with a constant birth rate during 10 Myrs, and 
detect 16602 neutron stars as the radio pulsars with a characteristic age of r < 5Myr. 
Scaling the simulated population to the observation, ~ 800 radio pulsars, provides ~1.3 
per century as the predicted birth rate. Figure [2] compares the distributions of the various 
characteristics of the simulated (unshaded histograms) and observed pulsars (shaded 
histograms) with a characteristic age of r < 5 Myr. We can see in Figure |2] that the present 
simulation reproduces well the observed distributions for the period time derivative, the 
spin down age and the surface magnetic field, although the simulation predicts more pulsars 
close to the Earth and more pulsars with a smaller radio flux, as bottom two panels show. 
The difference between the results of the observation and of the simulation may be caused 
by the assumptions within the present simulation. On the other hand, we expect that 
an adjustment of the slight difference between the simulated and observed population in 
Figure |2] will not change much the results of population on the 7-ray pulsars discussed in 
the following sections. 

3.1. Population "bright" 7-ray pulsars 

First, we sample only bright 7-ray pulsars, because the 7-ray pulsars with a larger 
flux measured on the Earth will be preferentially discovered by the Fermi 7-ray telescope. 
Moreover, because the detection of the bright point sources might not be affected much by 
7-ray background emissions, it might be true that most of (or all) bright 7-ray pulsars have 
already discovered by the Fermi, implying the distribution of the bright 7-ray sources will 
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be able to used to test the emission model. In this paper, we define "bright" 7-ray pulsars 
with a flux larger than F 7 > 10~ 10 erg/cm 2 s. In the first Fermi pulsar catalog (Abdo et al. 
2010a), the number of the bright 7-ray pulsars is 12 for the radio selected and 13 for the 
7-ray-selected pulsars (including the Gemniga pulsar). Later, Saz Parkinson et al. (2010) 
reported new eight 7-ray pulsars discovered in blind frequency searches. Four out of eight 
pulsars show the 7-ray emissions brighter than F 1 > 10~ 10 erg/cm 2 s. 

In Table 2, we show the simulated population of the radio-selected N rF>w -io and 
7-ray-selected N giF> i -w bright 7-ray pulsars. We assumed the beaming factor of the 7-ray 
emissions are / 7 = 0.4 and the birthrate is 1.3 per century. We present the results for K = 1 
in equation (|29|) . the maximum gap thickness of f max = f cr u and of f ma x = 1- 

As Table 2 shows, we predict that there are N rF>w -io = 18 ~ 23 radio-selected 
and N g)F>10 -w = 26 ~ 34 7-ray-selected bright 7-ray pulsars. Although the ratio of the 
simulated 7-ray-selected to radio-selected bright sources, 1.3-1.5, is consistent with the 
present Fermi bright pulsars, 17/12 ~ 1.4, our simulation predicts more number than 
the observed number (12 radio-selected and 17 7-ray selected pulsars). To explain this 
difference, we may expect that although all bright sources are actually discovered by the 
Fermi, the detection of the pulsations were failed due to the technical problem, which may 
be caused by the strong background etc. The recent Fermi first source catalog (Abdo et al. 
2010d_|) includes about 50 unidentified sources with a flux F 1 > 10 -10 erg/cm 2 s, and most 
of them can be categorized as a steady source. The present model indicates that 15 ~ 28 of 
the unidentified Fermi bright sources originate from the canonical 7-ray pulsars. 

Figure [3] compares cumulative distributions of various characteristics for the simulated 
and observed bright 7-ray pulsars, including the radio-selected and 7-ray-selected 7-ray 



1 http : //fermi. gsfc.nasa.gov/ssc/data/access/lat/lyr_catalog 
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pulsars. The results of the simulation are for f max = f crit . We discuss the combined 
cumulative distribution of the radio-selected and 7-ray-selected pulsars, (1) because we 
increase the number of observation sample (12+17=29) to increase the statistical accuracy, 
and (2) because some 7-ray-selected pulsars have been identified in radio bands (Camilo et 
al. 2009). 

We performed a Kolmogorov-Smirnov (KS) test to compare the two cumulative 
distributions for the bright sources. In Figure |3l we present the maximum deviation (D max ) 
between the two distributions and the p-value (Pk s ) of the KS-test. It can be seen that the 
hypothesis that the two distributions are drawn from the same parent population can not 
be rejected at better than ~ 80 % confidence level for the period, spin-down age and the 
surface magnetic field, than ~ 50 % for the period time derivative and the flux, and than 
~ 20 % for the distance, indicating those model distributions will be consistent with the 
observations. 



3.2. Comparison with Fermi 7-ray pulsars 

In this section, we present simulated population of the 7-ray pulsars applying the 
Fermi sensitivities of the six-month observations, and compare the simulated population 
with all known canonical 7-ray pulsars listed in Abdo et al. (2010a) and Saz Parkinson et 
al. (2010). For the radio-selected 7-ray pulsars, the Fermi sensitivities of the six- month 
observation is approximately described by F = 2 x 10 -8 photons/cm 2 s for the Galactic 
latitude |6| > 5° and 6 x 10~ 8 photons/cm 2 s for |6| < 5°, while for the 7-ray-selected pulsars, 
F = 4 x 10~ 8 photons/cm 2 s for |6| > 5° and 1.2 x 10~ 7 photons/cm 2 s for \b\ < 5°. 

Figures H] and |5] show the distributions for the various characteristics for the 
radio-selected and 7-ray-selected 7-ray pulsars, respectively. The shaded and un-shaded 
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histogram show the results for the Fermi pulsars and simulated pulsars, respectively. The 
maximum deviation of the cumulative distributions (D max ) and the p-value of the K-S 
test (Pks) ar e also indicated in each panel. The results are for f max = f crit . We can see 
that except for the distributions of the period time derivative of the radio-selected 7-ray 
pulsars, p-vale of the K-S test is larger than Pk s ~ 0.15, indicating hypothesis that the two 
distributions are drawn from the same parent population can not be rejected at better than 
~ 85 %. Therefore, we would conclude that distributions of the our simulation are not 
conflict with the Fermi observations. 

Although KS-statistic indicates the consistency between the simulated and the observed 
distributions, one can see an excess in the observed distribution (shaded histograms in 
Figures H] and |5]) compared with the simulated distribution at rotation period P ~ 0.1 s, the 
spin down age r ~ 10 5 yr and the distance d ~ 1 kpc. This may be cause by intrinsic feature 
of the pulsars or selection effect or enhancement of the birth rate due to the local Galactic 
structure. Because the excess is caused by only three or four pulsars, all possibilities can be 
considered. If the excess of the pulsars is due to the local structure, the Gould Belt may be 
possible site to produce the pulsars which make the excess in the distributions. For example, 
Cheng et al. (2004) discussed that about 20 7-ray pulsars with a flux F 1 > 7 x 10~ n are 
associated with the Gould Belt. The Gould Belt has an ellipsoidal ring with semi-major 
and semi-minor axes equal to ~0.5 and 0.35 kpc, respectively (Guillout et al. 1998), 
although a smaller size is also predicted (c.f. Perrot & Grenier 2002). The distance to the 
center of the Gould Belt from the Sun is ~ O.lkpc, indicating the neutron stars are born 
about 0.25-0.6 kpc away from the Sun. If we use the typical velocity of the neutron star 
~300 km/s, the displacement during r ~ 10 5 yr is d ~ 300 km/s x 10 5 yr ~ O.lkpc. Such a 
small displacement will not be able to explain the excess of the distance at 1 kpc, although 
the size of the Gould Belt can play an important role. We expect that if the excess in the 
distributions are caused by the selection effects or by the effect of Gould Belt, the excess 
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will decrease as increase of the 7-ray pulsar population. On the other hand, if it is the 
intrinsic feature of the pulsar population, the excess will remain even the population of the 
pulsars increase. 

In Table 2, we show the number of the simulated radio-selected (N rto t) and 7-ray- 
selected (N gtot ) 7-ray pulsars. Our simulation predicts that with the present sensitivity, the 
Fermi can detect N rtot = 40 — 61 of the radio-selected 7-ray pulsars and N gto t = 36 — 75 
of the 7-ray-selected 7-ray pulsars, which are more than the present Fermi pulsars, that 
is, 21 radio-selected and 24 7-ray selected pulsars (Abdo et al. 2010a; Saz Parkinson et 
al. 2010). To explain the difference between the observation and the model, the several 
reasons are expected; because most of the canonical 7-ray pulsars are located in the 
Galactic disk (c.f. Figure [TO]) . 1) the simulated sensitivity of the Fermi is too simple, 
namely, the sensitivity becomes much worse at the Galactic plane (\b\ ~ 0), 2) a strong 
7-ray background prevents the detection of the pulsed emissions from some pulsars, and (3) 
the 7-ray emissions from the pulsars may be missed by source confusion with the complex 
regions (e.g. Cygnus and Carina that lie on tangents to spiral arms), and/or with the 
unresolved sources that are not modeled in the diffuse backgrounds. 

As Table 2 shows, our simulation also predicts that if the Fermi sensitivity increase 
factor of two from the present one, Fermi can in principle detect N rto t x2 = 58 — 89 
radio-selected and N gtot X 2 = 79 — 152 7-ray-selected 7-ray pulsars. 

3.3. Implication for future observations 



Finally, we discuss model implication on future observations. In the observation, 
sensitivity of the detection of 7-ray emitting pulsars really depend on the observed direction 
(that is, high or Galactic latitude) and 7-ray-selected or radio-selected pulsars, implying the 
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present observed distribution of the pulsar characteristics may not represent the intrinsic 
population of the 7-ray pulsars that irradiate the Earth. 

Figure [6] summarizes number of the pulsars detected within the simulation as a 
function of the 7-ray threshold energy flux. Here, because we do not take into account 
the dependency of the sensitivity on the observed direction and on the radio-selected or 
7-ray-selected pulsars, Figure |6] presents the intrinsic number of the pulsars irradiating 
the Earth with a flux larger than the threshold flux. The solid and dashed lines represent 
numbers for the radio-selected and 7-ray-selected 7-ray pulsars, respectively. In addition, 
the thin and thick lines are results for the fractional gap thicknesses of f max — fcrit an d 
fmax = 1) respectively. In the simulation, for example, 18-23 radio-selected (solid lines) 
and 27-35 7-ray-selected (dashed line) 7-ray pulsars irradiate the Earth with a flux large 
than F 7 = 10~ 10 erg/cm 2 s measured on the Earth, while 64-89 radio-selected and 340-441 
7-ray-selected 7-ray pulsars irradiate the Earth with a flux larger than F 7 = 10 -11 erg/cm 2 s. 
Our simulation therefore implies that although the minimum energy flux of the known 
7-ray pulsars is approximately F 1 ~ 10~ n erg/cm 2 s, the Fermi has missed most of the 
7-ray pulsar with a 7-ray flux larger than F 7 = 10~ n erg/cm 2 s. 

Figures [7] and [8] show the distributions (unshaded histograms) of characteristics of the 
simulated radio-selected and 7-ray-selected 7-ray pulsars, respectively, with a 7-ray flux 
larger than F 7 > 10 -11 erg/cm 2 s. In the figures we also show the observed distributions 
(dashed histograms) by the Fermi. For the rotation period, the simulation predicts that 
the future observations will detect the radio-laud pulsars with longer rotation period 
(P > 0.5 s) as Figure [7] shows, if the sensitivity is improved. For the 7-ray-selected pulsars, 
on the other hand, the future observations will detect more pulsars with a shorter rotation 
period than 0.1 s. 

We find in Figures [7] and [8] that most of the simulated pulsars distribute with a distance 
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around d ~ 10 kpc and a flux F 7 ~ 10~ n erg/cm 2 s, while the distributions of the Fermi 
pulsars have a peak at d ~ 1 kpc and F 7 ~ 10 -10 erg/cm 2 s. Therefore, our simulation 
predicts that as the sensitivity is improved, more and more 7-ray pulsars with a distance 
lager than d ~ 1 kpc and a flux smaller than F 1 ~ 10~ 10 erg/cm 2 s will be detected by the 
future observations. 

We note that as Figures [7] and [8] show, our simulation indicates the peak of distributions 
for the period and the period time derivative come to P ~ 0.1 — 0.3 s and P x 5 x 10 -14 s/s, 
respectively, indicating the 7-ray pulsars with a rotation period of P ~ 0.1 ~ 0.3 s 
and a period time derivative of P ~ 5 x 10~ 14 s/s may be preferentially detected from 
the unidentified Galactic sources. This information might be useful to narrow down 
the searching range of the pulse period from the unidentified sources to detect new 
7-ray-selected pulsars. 

4. Discussions 

4.1. Ratio between radio-selected and 7-ray-selected pulsars 

In this section, we discuss the ratio between the radio-selected and 7-ray-selected 7-ray 
pulsars. With Table 3, we demonstrate how the ratio depend on the radio surveys. The first 
column shows the radio surveys used in the simulation, where M2, A2, and P2 represent for 
Molonglo 2, Arecibo 2 and Parks 2 surveys. "All" presents results with all radio surveys 
listed in Table 1. In addition, we increase factor of two for the sensitivity of each survey. 
The second and third column show the number of the radio-selected (N r ) and the 7-ray 
selected pulsars that irradiate the Earth with a flux larger than F 7 = 10~ 10 erg/cm 2 s 
and 10 -11 erg/cm 2 s, respectively. The fourth column shows the results using the Fermi 
sensitivity of six-month observations. We present the results for f max = f crit . 
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As Table 3 shows the ratio between the radio-selected and radio-quiet 7-ray pulsars 
depend on the sensitivity of the radio survey. For example, the present simulation predicts 
the ratio N g /N r = 4 using only Arecibo 2 and Parks 2 surveys, which will be consistent with 
the Ng/N r = 37/11 ~ 3.4 predicted by Cheng & Zhang (1998). Because radio emissions 
from some 7-ray selected pulsars can be detected by more sensitive radio survey, the ratio 
using "All" radio survey listed in Table 1 becomes N g /N r ~ 1.4, which will be consistent 
with the present Fermi result N g /N r ~ 17/12 ~ 1.4. 

Table 3 indicates that the ratio depends on the threshold energy flux, that is, the ratio 
increases as the threshold energy flux decreases, as second and third columns show. This 
tendency depends on the sensitivity of the radio survey. For example, because the "bright" 
7-ray sources are located relatively near the Earth, the radio emissions from those pulsars 
can be observed by the present radio surveys if the radio emissions point toward the Earth. 
On the other hand, dimmer 7-ray sources (F 1 ~ 10 -11 erg/cm 2 s) tend to have a more 
distance as indicated by bottom panels of Figures [7] and El The flux of radio emissions 
from those distant 7-ray pulsars may be under the sensitivity of the radio surveys, and the 
pulsars are categorized as 7-ray-selected pulsars. Therefore, the ratio of the 7-ray-selected 
and radio-selected pulsars increases as the threshold energy flux decreases. In fact, the 
present simulation predict that if we could detect all pulsars that irradiate the Earth with 
the radio emissions, the ratio of the 7-ray-selected and radio-selected pulsars does not 
depend on the threshold energy flux. 

Applying the sensitivity of the Fermi (fourth column), because the sensitivity for 
the radio-selected 7-ray pulsars is better than the 7-ray-selected pulsars, the ratio N g /N r 
is in general smaller than those of second and third columns, whose sensitivities for the 
radio-selected and 7-ray-selected 7-ray pulsars are the same each other. With the all radio 
surveys, the ration N r /N g ~ 0.9 will be consistent with the six-month Fermi observations, 
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N r /N g ~ 24/21 ~ 1.1. 

4.2. Comparison with unidentified Fermi sources 

Finally, we discuss the distributions of the 7-ray pulsars in the galaxy. As we discussed 
in section 13. 3[ the present model predicts most of the canonical 7-ray pulsars with a flux 
smaller than F 7 < 10~ 10 erg/cm 2 s have been missed by the Fermi telescope. On the 
other hand, the Fermi has found several hundreds of the unidentified points sources, 
implying some of them must be the canonical 7-ray pulsars. It will be useful to compare 
the Galactic distributions for the simulated 7-ray pulsars and the unidentified sources 
to have a hint for the origin of the unidentified sources. Figures [9] and HOI respectively, 
show the Galactic longitude and latitude distributions for the simulated 7-ray pulsars with 
Fj > 10~ n erg/cm 2 s (solid line) and the unidentified Fermi sources (dashed line). We 
selected the unidentified sources with the criteria that the observed 7-ray flux is larger than 
10~ n erg/cm 2 s and the variability index (V) in the catalog smaller than 23.21, larger than 
which indicates less than a 1% chance of being a steady source (Abdo et al. 2010c). In 
addition, we chose the unidentified sources with the curvature index, which is defined in the 
first catalog, larger than 5. A curvature index larger than 11.34 indicates a less than 1 % 
chance that the power-law spectrum is a good fit for 100 MeV-lOOGeV. The first catalog 
show thats the 7-ray pulsars typically have a curvature index larger than 5. 

We find in Figure M that the Galactic longitude distribution of the simulated 7-ray 
pulsars is qualitatively consistent with that of the Fermi unidentified sources, that is, both 
distributions have a peak around I = 0° and become minimum around I ~ 180°. In fact, we 
can see that the distributions of the existing radio pulsars younger than the characteristic 
age of r < 5 Myr is also qualitatively consistent with that of the unidentified sources. 
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In Figure (TUJ on the other hand, we can see that most of the simulated 7-ray pulsars 
are located around the Galactic disk, that is, |6| < 5° ~ 10°, while significant fraction of 
the unidentified sources are located above the Galactic disk. Because we could see that 
the longitudinal distribution of the simulated 7-ray pulsars is consistent with that of the 
existing radio pulsars younger than the characteristic age r <5 Myr, our model prediction 
is that most of high-latitude unidentified sources will not associate with the canonical 
7-ray pulsars. If the majority of the unidentified sources is originated from the Galactic 
sources, the millisecond pulsars may be possible candidate. Because millisecond pulsars are 
in general older than the canonical 7-ray pulsars and have a slower proper motion than the 
canonical pulsars, they are bounded by the Galactic potential. It is expected that a more 
fraction of 7-ray emitting millisecond pulsars is located at the higher Galactic latitude. On 
these ground, it will be important to perform a population study of the millisecond 7-ray 
pulsars to reveal the origin of the Fermi unidentified sources. The results of the population 
study for the millisecond pulsars using a Mont-Carlo technique will be discussed in the 
subsequent papers. 



5. Conclusion 

Inspired by the drastically increase of the 7-ray emitting pulsars by the Fermi 
telescope, we performed a population study for the canonical 7-ray pulsars with a 
Monte-Carlo method. We applied the outer gap model with a switching of the gap closure 
process from the photon-photon pair-creation to the magnetic pair-creation model suggested 
by Takata et al. (2010). Our simulation predicts that there are 18-23 radio-selected and 
26-34 7-ray-selected 7-ray pulsars with a flux lager than 10~ 10 erg/cm 2 s measured on the 
Earth. Applying Fermi sensitivity of six-month observations, our simulation predicts that 
40 — 61 radio-selected and 36-75 7-ray-selected pulsars can be detected. We showed that the 
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simulated distributions for the various characteristics (e.g. rotation period, characteristic 
age) can successfully explain those of the Fermi 7-ray pulsars. Our simulation also predicts 
that there are ~ 64 radio-selected and ~ 340 7-ray-selected 7-ray pulsars irradiate the 
Earth with a 7-ray flux larger than F 7 > 10~ n erg/cm 2 s. This indicates that most of the 
7-ray pulsars have been missed by the present observations, although their fluxes are larger 
than the minimum 7-ray flux detected by the Fermi so far. In particular, our simulation 
predicts that if the sensitivity of the instrument is improved, more 7-ray pulsars located 
with a distance larger than 1 kpc and a flux F 7 ~ 10 -11 erg/cm 2 s will be detected by in the 
future observations. The ratio of the simulated 7-ray-selected to radio-selected 7-ray pulsars 
is ~ 1.4 for F 7 > 10~ 10 erg/cm 2 s, while ~ 1 using 6-month Fermi sensitivity, although the 
ratio depends on the sensitivities of the radio surveys. Finally, we compared the Galactic 
distributions between the simulated 7-ray pulsars and unidentified Fermi sources. We 
saw that although the Galactic longitude distribution of the simulated 7-ray pulsars is 
qualitatively consistent with that of the Fermi unidentified sources, the simulated latitude 
distribution can not explain the observations, in particular for the higher latitude sources. 
This may suggest that the population study for the 7-ray emitting millisecond pulsars may 
help to understand the origin of the high-latitude Fermi unidentified sources. 
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Parks MB 


0.735 


21 


21 


1374 


2100 


0.250 


285 


3 


[260,50] 


[-5,5] 


Manchester et al. (2001) 


Jordell Bank 2 


0.4 


6 


40 


1400 


524 


1 


10 


5 


[355,105] 


[-1,1] 


Clifton et al. (1992) 


Swinburne IL 


0.64 


15 


21 


1374 


265 


0.125 


288 


3 


[260,50 


[5,15] 


Edwards et al. (2001) 



Table 1 : Parameters of each radio survey. 
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Number N r>F>10 -io N gtF>10 -io 


N r tot Ngtot 


Nrtot x2 -^tot x2 


T , fmax 1 

K=l 


23 34 


61 75 


89 152 


fmax fcrit 


18 26 


40 36 


58 79 



Table 2: N r)F>w -io and iV 9>F>10 -io are the number of the simulated radio-selected and 7- 
ray-selected pulsars with a flux of F > 10 -10 erg/cm 2 s, respectively. N rtot and N gtot present 
number of the simulated radio-selected and 7-ray-selected pulsars using Fermi sensitivity 
of the six-month observations; for the radio-selected pulsars, F = 2 x 10~ 8 photons/cm 2 s 
for the |6| > 5° and 6 x 10~ 8 photons/cm 2 s for \b\ < 5°, while for the 7-ray-selected pulsars, 
F = 4 x 10~ 8 photons/cm 2 s for |6| > 5° and 1.2 x 10 -7 photons/cm 2 s for |6| < 5°. For 
N r totx2 and N gtotx2 , we increases factor of two for the Fermi sensitivity. 
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Survey 


N r 




N g /N r 


N r 




N g /N r 


N r 


N g N g /N r 














F 7 


> io~ 10 


erg/cm 2 s 


F 7 


> io- n 


erg/cm 2 s 


Fermi sensitivity 


A2, P2 


9 


35 


3.9 


24 


372 


15.5 


17 


48 2.8 


M2, A2, P2 


13 


28 


2.2 


36 


367 


10.2 


26 


40 1.5 


All 


18 


26 


1.4 


64 


340 


5.3 


40 


36 0.9 


All (x2) 


20 


26 


1.3 


93 


319 


4.2 


47 


34 0.72 



Table 3: Ratio between the radio-selected (N r ) and 7-ray-selected (N g ) pulsars vs. radio 
surveys. The first column is list of survey, the second and the third columns present the 
number of the pulsars detected within the simulation with F 1 > 10~ 10 erg/cm 2 s and F 7 > 
10 -11 erg/cm 2 s, respectively. The final column shows the number of the pulsars detected 
within the simulation using the Fermi sensitivity of the six-month observation. The results 
are for the beaming fraction of / 7 ~ 0.4 and the birth rate of 1.3 per century. M2=Molonglo 2, 
A2=Arecibo 2 and P2=Parks 2. "AH" represents all surveys listed in Table 1. In 'All(x2)", 
we increase factor of two for the sensitivity of each survey. 
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Fig. 1. — Radio flux threshold for each radio survey as a function of the rotation period. The 
results are for T s k y = 150 K at 408 MHz and the dispersion measure of DM = 200 cm _3 pc. 
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Fig. 2. — Distributions of various characteristics of the radio pulsars with a characteristic 
age (r = P/2P) younger than 5 Myr. The shaded histograms and unshaded histograms 
present results for the observation and simulation, respectively. The data were taken from 
Manchester et al. (2005). 
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Fig. 3. — Comparison of the cumulative distributions for the various characteristics of 
"bright" 7-ray pulsars, F 7 > 10~ 10 erg/cm 2 s, including both radio-selected and 7-ray- 
selected 7-ray pulsars. The solid and dashed lines represent the distributions for the obser- 
vations and for the simulation, respectively. The p-value {Prs) of the Kolmoglov-Smirov 
test and the maximum difference (D max ) between two distributions are also indicated in each 
panel. The data were taken from Abdo et al. (2010a) and Saz Parkinson et al. (2010). 
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Fig. 4. — The distribution of the various characteristics for the radio-selected 7-ray pulsars 
for six-month Fermi observations. The shaded and unshaded histograms presents results 
for the observations and for the simulation, respectively. In addition, D max and Pks show 
the maximum difference of the cumulative distributions and the p- value of the KS-test, 
respectively. In the simulation, the sensitivity of the 7-ray detection is set at F = 2 x 
10~ 8 photons/cm 2 s for |6| > 5° and 6 x 10~ 8 photons/cm 2 s for \b\ < 5°. The data were taken 
from Abdo et al. (2010a) and Saz Parkinson et al. (2010). 
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Fig. 5. — The same with Fig. HJ but results for the 7-ray-selected 7-ray pulsars. In the 
simulation, the sensitivity of the 7-ray detection is set at F = 4 x 10 -8 photons/cm 2 s for 
\b\ > 5° and 1.2 x 10~ 7 photons/cm 2 s for \b\ < 5°. 
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Fig. 6. — The number of radio-selected (solid) and 7-ray-selected (dashed) pulsars detected 
within simulation as a function of the threshold energy flux. The thick and thin lines are 
results for the maximum fractional gap thickness of f max = f cr u and f max = 1, respectively. 
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Fig. 7. — The simulated distribution (unshaded histograms) of the characteristics for the 
radio-selected 7-ray pulsars that irradiate the Earth with a flux F 7 > 10 -11 erg/cm 2 s. The 



shaded histograms present results for the Fermi observation. 
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Fig. 8. — The same with Fig. but results for the 7-ray-selected pulsars. 
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Fig. 9. — The Galactic longitude distributions of the simulated 7-ray pulsars (solid line) and 
the Fermi unidentified sources with F 7 > 10" 11 erg/cm 2 s. The data were taken from the 



Fermi first catalog (Abdo et al., 2010c) with the criteria that the variable index less than 



23.11 and the curvature index larger than 5. 
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Fig. 10. — The Galactic latitude distributions of the simulated 7-ray pulsars (solid line) and 
the Fermi unidentified sources with F 7 > 10~ n erg/cm 2 s. The data were taken from Abdo 



et al. (2010c). 



